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A new family of j4jv-type Dunkl operators preserving a polynomial subspace of finite dimension 
is constructed. Using a general quadratic combination of these operators and the usual Dunkl oper- 
ators, several new families of exactly and quasi-exactly solvable quantum spin Calogero-Sutherland 
models are obtained. These include, in particular, three families of quasi-exactly solvable elliptic 
spin Hamiltonians. 



I. INTRODUCTION 

b . 

■ In the early seventies, Calogero Q and Sutherland [^6[ ^ introduced the celebrated exactly 

solvable (ES) and integrable quantum many-body problems in one dimension that bear their names. 
• These papers had a profound impact in the whole physics community, as reflected by the vast 

1^ I amount of literature devoted to the study of the mathematical properties and applications of 

these models. Among the most recent ones we could mention soliton theory p9, k2[, orthogonal 



' polynomials [fl, [l3, |27t , fractional statistics and anyons ||], random matrix theory [p9|, and Yang- 



> 



Mills theories~[|liy 20|| , to name only a few. Later on, Olshanetsky and Perelomov [^9[ explained the 
integrability of the original Calogero-Sutherland (CS) models by relating them to the root system 
of Aat type. These authors then constructed new families of integrable many-body Hamiltonians 
associated with all the other root systems. Furthermore, they showed that the most general 
I interaction potential for these models is proportional to the Weierstrass p function. 

. A considerable effort has been devoted over the last decade to the extension of CS models to 

' particles with spin. These models are a step forward towards the unification of the CS scalar models 

and integrable spin chains, like the Haldane-Shastry model [|2^, Several different techniques 
' have been used to construct spin counterparts of the scalar CS models, including the exchange 

I operator method |]3l| , the Dunkl operators formalism ^ , the supersymmetric approach , and 

reduction by discrete symmetries [ p3| . In the quantum case, only the rational and trigonometric 
pi il (or hyperbolic) spin CS models have been constructed, both in their An |[ ^, 31, 



D . and BCjy |Q versions. In the An case, the integrability and the exact-solvability of these models 

' both follow from the fact that the Hamiltonian is related to a quadratic combination of some family 

of Dunkl operators. 
■ The Dunkl operators 

T, = — +aJ2^—i'^-K^J), 1^1,. ..,N, (1) 

were originally introduced in in connection with the theory of orthogonal polynomials associ- 
ated with finite reflection groups. In the latter expression, a is an arbitrary real parameter and 
the sum runs over j = 1, . . . ,i — l,i + I, . . . , N. The permutation operators Kij — Kji act on an 
arbitrary function /(z), with z = (zi, . . . , zn) S K"'^, as 

{Kijf) (zi, . . . , Zi, . . . , Zj, . . . , zat) = /(zi, . . . , Zj, . . . , Zi, . . . , zn). (2) 

Using the relations 

= 1, KijKjk = K^kKij = KjkKik, K^jKki = KkiKij, (3) 

where k, I take different values in the range 1, . . . , A^, one can establish the commutativity of 
the operators (|lj) and prove that Ti,Kjk, i,j,k = 1,...,A^, span a realization of a degenerate 
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Hecke algebra (see |^ for more details). Since the rational spin CS Hamiltonian is related to a 
polynomial in the Dunkl operators (|l|) , these operators yield a complete set of commuting integrals 
of motion. In addition, the spectrum of the Hamiltonian follows immediately from that of the 
Dunkl operators, which can be easily computed ||. 
The previous considerations also apply to the operators 

= A + a V (1 - K,,) + a V (1 - K,,) + 1 - z, (4) 

OZi Zi — Zj Zi — Zi 

j<i ■> j>i ■> 

i — 1, . . . , N, introduced by Cherednik ^ in connection with the trigonometric spin CS model. In 
other words, the operators Ti commute, have an easily computable spectrum, and can be used to 
obtain a complete set of integrals of motion and the spectrum of the Hamiltonian. It has become 
customary in the literature to refer to both families of operators Ti and T; as Dunkl operators. 
Recently, some partiallysolvable deformations of the scalar CS models with an external potential 



have been proposed 24, For these models — in contrast to the CS models listed in [g9| — 
only a finite-dimensional subset of the spectrum can be computed algebraically. Following Turbiner 
and Ushveridze [^0| we shall use the term quasi- exactly solvable (QES) to refer to this type of 
models; see also the reviews In all these models, the Hamiltonian can be expressed as 

a quadratic combination of the generators of a realization of 5l{N + 1) by first-order differential 
operators preserving a finite-dimensional space of smooth functions. The action of the Hamiltonian 
in this space can thus be represented by a finite-dimensional constant matrix. The eigenvalues of 
this matrix belong to the spectrum of the Hamiltonian provided the corresponding eigenfunctions 
satisfy some appropriate boundary conditions. The Lie algebra s[(7V -I- 1) is usually referred to as 
a hidden symmetry algebra of the Hamiltonian in these models. 

In this paper, we propose a general procedure for constructing (Q)ES spin CS models, close in 
spirit to the hidden symmetry algebra approach to scalar QES models. The starting point of our 
construction is the well-known fact that the two standard families of Dunkl operators (|l|) and 
admit an infinite sequence of invariant polynomial subspaces of finite dimension. One of the main 
novelties in our approach consists in the introduction of a new family of commuting Dunkl operators 
which, together with the other two families (|l]) and (^), is shown to preserve a single polynomial 
subspace of finite dimension. We then prove that certain quadratic combinations involving all 
three families of Dunkl operators always yield a spin CS Hamiltonian. The QES character of these 
Hamiltonians follows immediately from the fact that the Dunkl operators admit a finite-dimensional 
invariant subspace. Moreover, if the original quadratic combination does not involve the new family 
of Dunkl operators, the resulting Hamiltonian preserves an infinite sequence of finite-dimensional 
subspaces of smooth functions, which we shall take as the definition of exact solvability. The 
linear space spanned by all three types of Dunkl operators is then shown to be invariant under 
the projective action of the group GL(2,R). We make use of this fact to perform a complete 
classification of the resulting (Q)ES spin CS models. All the previously known exactly-solvable 
spin CS models of A^^ type appear as particular cases, arising from a quadratic combination of 
a single type (either (|l|) or (Q)) of Dunkl operators. In addition, we obtain many new spin CS 
models, both exactly and quasi-exactly solvable. These include, in particular, several elliptic QES 
spin CS models. To the best of our knowledge, these are the first examples of solvable quantum 
spin CS models involving elliptic functions. 



II. A NEW FAMILY OF Ajv-TYPE DUNKL OPERATORS 

In this section we shall define a third family of Dunkl operators that preserve certain finite- 
dimensional polynomial subspaces. These three families shall be used in the following sections to 
construct exactly and quasi-exactly solvable spin many-body Hamiltonians. 

Let us begin by introducing the polynomial subspaces 

7?.,„(z) = span J]^ z'' : k < m, i = \, . . . ,n \ , (5) 



{N N \ 

1=1 i=l J 



(6) 
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which shall be referred to as the rectangular and triangular modules, respectively, by analogy with 
the two particle case |l5| . 

A well-known property of the Dunkl operators (|l|) , (^) is the fact that they preserve the triangular 
module 7'„(z) for arbitrary n. Seemingly less known, but central to our construction, is the fact 
that they preserve the rectangular module 7?.„i(z) for arbitrary m as well. On the other hand, the 
differential parts of the Dunkl operators (|l|), (jj) together with the differential operator zfdz^, span 
a realization of s[(2). Inspired by this fact, it is natural to suggest the following ansatz for a third 
set of Dunkl operators: 

d 

= ^"^Ti + /y (z)(l - Kij), i = l,...,N, 

ozi 

where m is an arbitrary non-negative integer and fij{z) is a function anti-symmetric in This 
new family does not preserve the module 7-'n(z), but for a suitable choice of the functions fij{z) it 
will be shown to preserve the module TZmi^). 
To this end, let us define the operators 



Q°-iE7^(l-^^.), (7) 



Zi Zj 
Zn 



where a is a real parameter and the sum runs over j = 1, . . . z — 1, i -I- 1, . . . iV. 
The following lemma will be important in the sequel: 

Lemma 1. For any non-negative integer n, the rectangular module 7?.„(z) is invariant under the 
action of the operators , and Qf . The triangular module Vn (z) is invariant only under the 
action of the operators and Qf . 

Proof. It suffices to prove that the inclusions 

^(l-i^„■)7^„(z)c7^„(z), e = ±,0. 



-3^(l-if„)7'„(z) C7'„(z), e = -,0, 

Zi Zj 

hold for any pair of indices i ^ where 

= 1 ' ^Ij = ^» + , h+ = ZiZj . 
The action of these operators on an arbitrary monomial 

hh . z^z^ - i'z"; ^ ; 



2» Zj fc^l Z^ h Z^'Z^ 

yields a polynomial, since z^^ z^ — z^ z^^ is a multiple of Zi — Zj. The homogeneous degree of this 
polynomial is either (if li — Ij) ox t ^ (i^ 7^ ^j)- Therefore, if the original monomial 

belongs to 'P„(z) the resulting polynomial is also in 7'„(z) for e = — ,0, but lies outside 7'„(z) for 
e +. 

On the other hand, the degrees of the variables Zk in the resulting polynomial remain equal to 
Ik if k ^ i, j , while the degrees of Zi and Zj satisfy 

deg(zfc) < niax(/i, Ij) -1 + d^, k ^ ij , 

where 
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Therefore, if the original monomial belongs to the space 7?.„(z) so does the resulting polynomial 
in all three cases e = ±, 0. □ □ 

The following three sets of Dunkl operators shall be the building blocks for the construction of 
several new (quasi-)exactly solvable spin CS models: 



where i = 1, . . . , iV, a is a real parameter, and m is a non-negative integer. Note that the operators 
J~ coincide exactly with the Dunkl operators (|l|) , while the operators differ from the Cherednik 
operators (Q) by a linear combination with constant coefficients of the permutation operators Kij , 
namely 

f.^J^ + lY.i^-K,,)^^Y.i^-K,,) + ^ + l-^■ (9) 

j<i j>i 

The operators J^' are, to the best of our knowledge, new. 

The operators and Kij obey the following commutation relations 

[x.„j^] = o, k,,j^ = j;k,,, (11) 

where e = ±, and the indices k are all different. The set of operators 

{J!,K,, : i,j = 1,...,N}, e = ±,0, 

spans a degenerate affine Hecke algebra, see This is clear for e = ±, while for e = 0, it follows 
from (||) and the commutativity of the Cherednik operators (^. 

The key property in our construction of (quasi-)exactly solvable spin CS models is the fact that 
the operators (|8[) possess invariant polynomial subspaces. 

Theorem 1. The operators and Jf preserve the modules Vn{z) and 7?.„(z) for an arbitrary 
non-negative integer n. The operators preserve the module TZm{z), but do not preserve the 
modules 7^„(z) and 7?.fc(z) for k ^ m. 

Proof. The statement follows from Lemma |l| and the fact that the differential parts of and J° 
preserve the modules 7^„(z) and TZ„{z) for any non-negative integer n, whereas the differential part 
of preserves the module TZkiz) only for k = m. □ □ 

The following corollary is an immediate consequence of Theorem |^: 

Corollary 1. Any polynomial in the operators leaves invariant the rectangular module TZ^^z). 
In addition, if the polynomial does not depend on , it preserves the modules 7?.„(z) and Vniz) 
for all n. 



III. CONSTRUCTION OF SPIN CALOGERO-SUTHERLAND MODELS 

In the previous section we have introduced a new set of Dunkl operators preserving the space 
of polynomials TZ,n- Here we shall make use of all three sets of Dunkl operators (^) to construct 
some multi-parameter families of spin CS models. 



5 



Consider the spin permutation operators 5*^, i,j = l,...,iV, whose action on a spin state 
|si, . . . , sjv), -M < s,< M, with M e iN, is given by 

Sij\si, ...,Si,...,Sj,...,SN)^ |si, ...,Sj,...,Si,...,SN). (12) 

Note that the operators Sij obey the identities ^ with Kij replaced by Sij. Let S denote the 
hnear space span { \si, . . . , sn)} The action of the operators Sij in 6 is thus represented 

by (2M + l)^-dimensional symmetric matrices. 

The starting point of our procedure is the following quadratic combination of the Dunkl opera- 
tors d): 

(13) 

where Cee',Ce, e, e' = ±,0, are arbitrary real constants. The term ^^i{Jj~ , Ji'} differs from 
X^iGZi)^ by a constant operator (see Appendix), and for this reason it has not been included 
in (p^. We emphasize that only the particular cases of ( p^ ) 



-^2 



have been previously discussed in the literature in connection with CS models; see ^, 13, 3^, Q 
and references therein. 

As it is customary, we shall identify Kij, Sij, and H* with their natural extensions Kij I, 
i(S)Sij, and H* (g) I to the tensor product C[zi, . . . , zn]^&. The following lemma is an immediate 
consequence of Eq. (|ll|) . 

Lemma 2. The (Q)ES differential- difference operator H* commutes with Kij and Sij for all 
i,j = 1,...,N. 

This property plays a crucial role in the construction of spin CS models; see, for instance, Ref. 

Let A be the projection operator on states antisymmetric under the simultaneous interchange of 
any two particles' coordinates and spins. In terms of the total permutation operators n,j = Kij Sij, 
the operator A can be alternatively defined by the relations A = —A, j > i = 1, . . . , N. Since 
Kfj = 1, these relations are equivalent to 

i^yA = -5,yA, j>i^l,...,N. (14) 

For the lowest values of N the antisymmetrizer A is given by 

= 2 : A = 1 - , 

= 3 : A = 1 - - - Has + nisHis + niaHsa . 

In general, A is an {N — l)-th degree polynomial in the total permutation operators Uij. It thus 
follows from Lemma ^ that H* commutes with A. 

Suppose that /(z) is an eigenfunction of H* with eigenvalue A. For instance, / could be one 
of the polynomial eigenfunctions that H* is guaranteed to possess in TZrn- Given any (constant) 
spin state |ct) G 6, the spin function ip — A[f{z)\a)] is also an eigenfunction of H* with the same 
eigenvalue A. 

Next, we introduce the matrix differential operator H obtained from H* by the formal substi- 
tutions Kij — !■ —Sij, i,j = l,...,N. The relations ([l^ ) imply that is a spin eigenfunction of H 
with eigenvalue A. 



Using the formulae ( Al)-([A6|) for the sums of the squares and the anticommutators of the Dunkl 



operators (jg) given in the Appendix, we obtain the following explicit expression for H: 
~H^J2 {Pi^^)9l + Q{zi)d,^ + R{z^)) + 2ac++(l - ™) ^ 

i i<j 



2 ^ 

■ {c++{zi + ZjY + ca+{zi + Zj) + cqo) Si-j + — cqoE ^ ^ij^ik), 

i<j i,j,k 
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where 



P{z) = C++Z* + co+z^ + cooz^ + co-z + c — , 



Q{z) = Q{z)+[l--)P'(z), 

Q{z) = c+z'^ + cqz + , b=l + m + a(N -1), 



(16) 



R(z) c++(6 + m(m - 2) - 1)2^ + 



6 + m(m — 1) — 1 



+ ^(2(6-l) + m(m-2))-^co 



and ^ ' denotes summation in i,j, k with i ^ j ^ k ^ i. 

In the final step, one performs a gauge transformation with a suitable scalar function /i(z), 
followed by a change of variables z = C(^)7 ^ — i^ii ■ ■ ■ 1 xn), 



(17) 



(18) 



z=C(x) 

in order to reduce the gauge spin Hamiltonian H to the Schrodinger form 

^ = -EC + nx), 

i 

where ^(x) is a Hermitian matrix-valued function. 

Proposition 1. The operator H in Eq. ( |l5| ) can be reduced to a matrix Schrodinger operator ( |l8| 
by a change of variables z = C(x) and conjugation by a scalar gauge factor fi{z). 

Proof. The gauge transformation ( [l7| ) with gauge factor 

N 

f,{z)^l[{z,-z,ri[p{. 

together with the change of variables 



exp 



^ Q{y) 



(19) 



dy 



z = l,. 



map the gauge spin Hamiltonian _ff to a matrix Schrodinger operator dlq), with potential 



(20) 



(Zi - Zj)(Zi - Zfe) 



Q(^.) 



l<j 



P{z,)+P{zj) 

{Zi - ZjY 



{a + S^j) 



2ac++(l - m-) ^ ZjZj - aE ('^++(^'' + ^if + c+o(2i + Zj) + coo) S"*: 



-E 

4i?(z.) 



^P'(z,)' + Qiz^f - 2Q{z,)P'{z,) ] ~ P"{z,) + 2Q\zi) 
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Coo 



E' (1 - 



□ 



z=C(x) 



(21) 



□ 



Remark 1. The gauge factor //(z) in (|l^) was introduced in [ p^ in connection with a generalization 
of the theory of QES models with one degree of freedom to many-body problems. The existence of 
a (matrix or scalar) gauge factor and change of coordinates reducing a given matrix second-order 
differential operator in iV > 1 variables to a matrix Schrodinger operator ( p^ is not guaranteed a 
priori. In the scalar case this problem was first addressed by Cotton pO[ |, while the matrix case 
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has been studied recently in Ref. . In fact, the quadratic combination H* has been chosen so 
that a scalar gauge factor and change of variables can be easily found for H. For instance, the 
term "^^IJ^ , J^] has been discarded because it involves first-order derivatives with matrix- valued 
coefficients, which are usually very difficult to gauge away. 

Remark 2. The change of variables (po|), and hence the potential V{x), are defined up to an 
arbitrary translation for each coordinate Xi, i = 1, . . . , N. We shall see that this arbitrariness can 
be removed in some cases by requiring the potential to be invariant under sign reversals of any 
coordinate Xi. 

If ip{z) is one of the eigcnfunctions of H with eigenvalue A constructed above, the spin function 
?/'(x) — /i(^(x))(p(^(x)) is clearly an eigenfunction of H with the same eigenvalue. Note, however, 
that we have not imposed so far any boundary conditions on the eigcnfunctions "0- In general, the 
parameters a, Cee', and Cg defining H* should satisfy certain constraints in order to ensure that the 
appropriate boundary conditions are satisfied. 

The following proposition is an immediate consequence of the previous considerations: 

Proposition 2. The spin Schrddinger operator (fsl) with potential (|2l|) leaves invariant the module 



= M(C(x))A(7^™(C(x)) ®6). (22) 
In addition, if C++ = co+ = c+ = 0, it preserves the modules Ain o-nd 

AA„ = a*(C(x))a(7'„(C(x))®6), (23) 

for any non-negative integer n. 
If we add a constant term 

Vo=jo+liY.S^j+-f2j2's^jS^k, 7»eR, (24) 

i<j i,j,k 

to the potential (^l]), the previous procedure for constructing eigcnfunctions of H still applies to 
H + Vo. Indeed, the associated operator {H + Vo)* is obtained by adding the term 

Vq* = 70 - 71 X! ^v^zk 

i<j i,j,k 

to the initial operator H*. Our assertion follows from the fact that Vq preserves the modules Vn 
and TZn for all n, commutes with all the permutation operators K^j^ and acts trivially on 6. This 
observation will be used in what follows to simplify the formula for T^(x) by dropping terms of the 
form (H). 



IV. CLASSIFICATION OF SPIN CALOGERO-SUTHERLAND MODELS 



We have seen in Section [II that any quadratic combination of the form (13) yields a spin CS 
model for which a number of eigenvalues and their corresponding eigcnfunctions can be computed 
in an algebraic fashion. In this section we shall obtain a complete classification of the potentials 
constructed in this way. 

The form of the potential V{x.) in Eq. ( pl|) depends on the choice of parameters Cee' and Cg, 
e, e' = ±, 0. The parameters Cee' which define the polynomial P are of particular significance, since 
they determine the form of the change of variables (|2^). However, different sets of parameters 
Cee', Ce defining the operator H* may give rise to the same potential. Indeed, there is a group 
of residual transformations preserving the vector spaces span{ J~, J°, J^^, i — l,...,iV. The 
image of the operator H* under these transformations is still of the form dl3), albeit with different 
coefficients c^^i and c^- We shall make use of this fact to classify the starting multi-parameter 
family of operators ( p^ ) into conjugacy classes. This will provide a complete classification of the 
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spin CS models obtainable within our framework. The ideas used for this classification are similar 
in spirit to those applied in Refs. [^8[ 19 to classify one-particle Lie-algebraic QES Schrodinger 
operators. 

Consider the mapping 

J|(w) ^ J^z) = J^(w(z)) /i-i(z) , (25) 

where w = {wi, . . . ,wn) is given by the projective action of GL(2,R) on RP^ (Mobius transfor- 
mation) 

z = l,...,7V, A = a5-/37^0, (26) 

and the gauge factor /im(z) is defined by 

N 

lim{z)=l[{^z, + 6y^. (27) 

4=1 

The following lemma is an immediate consequence of the definition of the Dunkl operators (^): 

Lemma 3. The mapping ( p5| ) acts linearly on the vector spaces span{ J^^, Jf, Ji~}: i = ■ ■ ■ i N , 
as 

'j;+(z)\ ^ fa' 2a/3 

J?(z) = ^ U7 aS + P^ (3S\ J0(z) I . (28) 
\j-{z)) 27<5 




It follows from the previous lemma that the operator H* defined by 

7?*=A^,„(z)i/*(w(z))M„i(z) 
is still a second degree polynomial in the Dunkl operators J|(z), whose coefficients c^e', can be 



easily computed using Eq. (|22). The corresponding gauge spin Hamiltonian H can be obtained 
from Eq. ( p^ ) by replacing the coefficients Cse' and Cg by their counterparts c^^i and c^. In particular, 
the polynomials P{z) and Q{z) in Eq. ( p!6| ) are replaced by the polynomials 

P{z) — c++z'^ + cq+z^ + cqoz'^ + cq-z + c — , Q{z) = c+z'^ + cqz + c_ . (29) 

Expressing the coefficients Cee' and Ce in terms of the original coefficients Cee' and Ce, one easily 
arrives to the explicit formulas 

Recall that the (irreducible) multiplier representation pn^i of GL(2, M) on the space of univariate 
polynomials of degree at most n is defined by the linear transformations 



p{z) ^ p{z) = A*(7z -f (5)"p 



az + (3 
7z + (5 



By Eq. (|3C|), the polynomials P and Q defining H transform according to the representations 
and p2,-ij respectively. Note also that the Dunkl operators (|^) transform according to 
the representation P2.-1; see Eq. (p8|). The (nonzero) orbits of the representation pi-2 can be 
parametrized by the following canonical forms ||l9|, IT] : 

±1, z, ±1/(1 -z^), ±1/(1 -Kz^), ±z/z^ ±z/(l + z2)2, 

±v{\^z'){i-ez'), ±v{i + z'){i + ez'), v{i^z'){i-e + pz'), 

where v > Q and Q < k < 1. The above list of canonical forms can be further reduced without any 
loss of generality using the complex projective (linear) transformation w — iz. Since the projective 
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transformation w = iz also induces the mapping Ce f—^ = i'^Cg for the coefficients of the polynomial 
Q, the resulting canonical form leads to a real potential provided that the initial coefhcicnts c± 
are purely imaginary. The reduced list of canonical forms will be conveniently taken as 

1) 1, 6) iy{l + z^)\ 

2) z, 7) z^){l-Pz^), 

3) iy{z^-l), 8) iy{z^ -l){l-k'^z^) , (31) 

4) 1/(1 -z2), 9) „{l-z^){k'^ + k^z^), 



2 



5) vz 

where /c'^ = 1 — fc^. By choosing P{z) in Eqs. (|l5|), ( p^ ) in each of the canonical forms (|3l|), one 
obtains a complete classification of the spin CS models with potential (pl|), which are (Q)ES by 
construction. Recall that in the one-particle case the canonical forms 1), 2), 3), 5) give rise to 
rational or hyperbolic potentials, while the remaining five yield periodic (trigonometric or elliptic) 
potentials jlH . 

Remark 3. The operator H* in (|l^) is easily seen to preserve the module 5m (z) of symmetric 
polynomials in zi, . . . , zat of degree at most to, on which the permutation operators Kij act as the 
identity. Hence the antisymmetrizer A acts on the space Sm ® & as the tensor product I (8) Aq, 
where Ag is the spin antisymmetrization operator, and the i7- invariant module fj,A(Sm ^ &) 
factors as the tensor product {^Sm) ^ (Aq©). The spin permutation operators Sij on this 
module reduce to —I. Therefore, the restriction of the Hamiltonian H to this space is simply 
(H\„ ) (g) I. Thus the scalar Schrodingcr operator L, leaves invariant the module 

fj,Sm- It follows that replacing the spin permutation operators Sij by —1 in one of the (Q)ES spin 
potentials listed below, one obtains a corresponding (Q)ES scalar potential. The scalar potentials 
so constructed include as particular cases all the potentials presented in ||l6|. 



For each of the canonical forms (|31|), the potential turns out to decompose as 

y(x) = ^C/(xO + l^nt(x), (32) 
i 

where U plays the role of a external field potential, and the interaction potential Vint is of the form 



Mnt(x) = {V'i^^ - ^j) + V^i^^ + Xj)) «(« + , (33) 



i<j 



with either = or = V~ . Indeed, making use of the identiti es (|A7| )-( A15| ) in the Appendix 
and recalling that co+ — for all the canonical forms, the potential (plj) reduces after some algebra 
to the form (|3^ ) with 

U = c++(l - b^)z^ + bc+z+ (3F'(z)2 + AQizf - 8Q(z)F'(z)) , 

Vint = 2^(Zi - z^y^ {c+j^zlz'j + CQQZiZj + ^ [zi + Zj) + c — ^ a{a + Sij) . 

i<j 

We have discarded here a constant term of the form (p^), in accordance with the observation 
at the end of Section The interaction potential takes the form ( ^3|) after performing the 
change of variables ( pO| ) and using some identities for the corresponding function Ci^)- Moreover, 
if the function is nonzero, it can be reduced to V~ by a suitable coordinate translation (see 
Remark 3). 

We shall now present the list of (Q)ES spin many-body potentials obtained from the canonical 
forms (p]). Note that the scaling (cee', Ce) t-^ (Acee', Ace) induces the mapping 

y(x;C<:<:',Ce) F(x; ACee', ACe) = AX^(VAx;Cee',C<:) 

/x(x;Cee',Ce) l-^ ^(x ; A C^' , A Ce) CX /L((VA X ; Qe' , Ce) 
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of the corresponding potentials and gauge factors. For this reason we shall list the potentials for a 
suitably chosen value of the parameter v in Cases 3-9, or a suitable multiple of P{z) in Cases 1,2. 
The notation 



shall be employed in what follows. 

Case 1. P(z) = -. 

4 

X 

Change of variables: = 2 • 
Gauge factor. 



External potential: 

U{x) = a]_x* + Aaoa+x^ + 4{al + 2a-a+) x^ + 2(8q;_q;o + 6a+) x . 
Interaction potential: 

14„t(x) = 2^(a;r.)-2a(a + 5,,). 

Case 2. P{z) = 4z. 
Change of variables: z = x^ . 
Gauge factor. 

i<j i 

External potential: 

U{x) = alx^ + 2aoa+x'^ + {al + a+{3b + 2a-))x^ + ^ ((2a_ - bf - l) - 
Interaction potential: 

M„t(x) = 2^] ((xr.)-^ + {x+y') a{a + Sij) . 



Case 3. P{z) = 4(^2 - 1). 
Change of variables: z = cosh 2a;. 
Gauge factor. 

At(x) = J|(sinha;^r sinha;+)'']J(sinha;i)5(i+"-^)(cosha;i)^(i+2«''-«-'') exp(^ cosh2a;i). 

i<j i 

External potential: 

U{x) = a\ cosh^ 2x + 2Q;+(ao + b) cosh2a; + 2(a+ + a_)(ao — b) cosh 2a; sinh~^2a; 
+ ((a+ + a-f + (ao - bf - l) sinh"^2a; . 

Interaction potential: 

yi„t(x) = 2 ^ (sinh-^a;-. + sinh-^a;+. j a{a + % ) . 



Case 4. P{z) = 4{1 - z"^). 
Change of variables: z = cos 2x. 
Gauge factor. 

/x(x) = nisina;,^ sina;+)"[|(sin:r,)5(i-"-'')(cosxi)5(i+«-2«o-*') exp( - ^ cos2a; 

i<j i 

External potential: 

U{x) = — Q!+ cos^ 2x + 2a_|_(6 — ao) cos 2a; + 2(a+ + a_)(6 + ao) cos 2a; sin~^2a; 
+ ((a+ + a-f + {b + aof - l) sin-22a; . 

Interaction potential: 



^int(x) = 2 ^ sin + sin \+-^ a{a + Sij) 

i<j 



Case 5. P{z) = Az^. 
Change of variables: z = e^^. 
Gauge factor. 

/x(x) = Yl sinh^Xj^ Yl ^ (a+e^^* - a_e~^^*^ + (ao - m)a 

External potential: 

U{x) = a\e^'= + 2a+(ao + 6) e^^ + 2a_ (ao - 6) e'^^ + a^.e""^ . 
Interaction potential: 

Mnt(x) = 2 ^ sinh~^a;r. a(a + Sij) . 

i<j 

Case 6. Piz) = {1 + z^f. 
Change of variables: z = tana;. 
Gauge factor. 

/j-(x) — JJsin°a;j^ JJcos™a;j exp (a+ + a_)a;i + - (a_ — a+) sin2xi — - aocos2xi 

i<j i 

External potential: 



U{x) = 2 ((°^+ ^ °^-)^ ^ Q^o) '^os4.7: + 2(Q:"i — a^ + 6ao) cos 2a; 



+ ao(a_ — a+) sin 4a; + 2(ao(a+ + a_) + 6(a+ — a_)) sin2x . 



Interaction potential: 



'(4nt(x) = 2 ^ sin "^Xij a(a + Sij) . 



Case 7. Piz) = A{1 - z^){l - Pz"^). 

Change of variables: z 
Gauge factor. 



dn2x 



Mx)=n(i 



sn x^^anx^^ 



k'^ sn'^x^J sii^x^j 



n 



(sn2a;j)3(i-''-:^)(dn2a;i)'^ 



X (dn 2a; j + cn 2a;i) ^"''^ (dn 2a;j + fc cn 2a;j) 
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Here (and also in Cases 8 and 9) the functions sna; = sn{x\k), cnx ^ cn(x|fc), and dnx = dn(.x|fc) 
are the usual Jacobian elliptic functions of modulus k, and k' = \Jl — k"^ is the complementary 
modulus. 

External potential: 

U(x) = A-r sn^2x + Br on 2x dn 2x + sn~'^2x (C7 + D-j on 2x dn 2x) , 



where 



A7 
Br 



2 

P2 



Dr=j^,{a^ + a.)(b+^) 



Interaction potential: 



yi„t(x)=2^ 

i<j 



+ 



a(a + S'jj) 



Case 8. P(^) = 4(^2 _ _ )t/2^2)_ 
Change of variables: z 



Gauge factor. 



dn2x 



(en 2a; j) 



.)5 [l-''-pi7(«++fc'c«o+fe"a_)] 



X sn 



2a;i) 5 (dn 2a;,)" (l + dn 2a;i) {k' + dn 2a;i) 



External potential: 



where 



U{x) = sn-22a; {Ag + Bg dn 2a;) + cn-22x (Cg + £>8 dn 2x) , 
^8 = ;^(a+ + a-)(|^-&), 

Cs = k'\b- - 1) + ^(if + 26) + ^ (a+ + k-a.f , 
D, = ^ia^ + k-a.)(^ + b). 



Interaction potential: 



Vint (X) =2^ 



dn^x. 



dn2-^+ 



2 - 2 - ' 2 + '2 + )a(a + 5y). 



i<3 ^ ^3 U U U 



Case 9. P{z) =A{l- z^){k'^ + k'^z'^). 
Change of variables: z = cn 2x. 
Gauge factor. 



i<3 



sn x,„. sn X J. dn x^^ dn x y -j-j- 
fc^sn^x"- sn^Xy / 



X (l + cn2xi)^^"++"-^exp 



(sn 2xi) ^(1-"-'') (dn 2xi) ^ (i+"o-&) 
(Acn2x,) 



k^a- - k'^a+ , / k 



2fcfc' 
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External potential: 



U{x) = (in ^2.T (Ag + Bg cn 2x) + sn (Cg + Dg cn 2x) . 



where 



Ag = fc'2(l - 62) + fc'2ao(2& - ao) + -1 (A:'2a+ - k^a_f , 

Bq ^ 2{b - ao) {k'^a+ - k^U-) , 
Cg = (6 + ao)^ + (a+ + a-)^ - 1 
Dg = 2(6 + ao)(a+ +a_) . 



Interaction potential: 

Vint (X) =2^ 



2 4 

cn a;/ 



sn^x,-,- dn^a;,^, sn^x,^ dn^x 



x— -ij -^i-j 

Remark 4- In Case 7, the alternative canonical form 



a{a + Sij) . 



P{z) ^ - g^z - 53 , 52 > 2753 . 

leads to a spin generalization of the QES potential involving Weierstrass functions studied in ||l^ . 
The corresponding change of variables is 2; = p{x + uj3), where p{x) = p(a;|52, 53) is the Weierstrass 
function with invariants 52, 93, and 2a;3 is its purely imaginary fundamental period. The gauge 
factor reads 



(p'(a;, + ^3))'^ 



X {p{xi + W3) - ei) ' (p(xi + ujz) - 62) {p{xi + ujz) - 63) 



\73 



where 



I3 



52a+ + 12(ejao + «-) 



(j, fc, Z) = cyclic permutation of (1, 2, 3) . 



24(ej - efe)(ej - e/) 

and Cj are the real (different) roots of P{z). The external and interaction potentials are given by 

U{x) = 4/3(/3 - l)p(2a;) + Ap{x + lo^) + {Bp^{x + CJ3) + Cp(x + W3) + £») (p'(a; + 1^3))"^ 
where 

4 

A = 2ao + -a+(2a+ + 36), 
9 

B = 4a2 + i (2a+ - 35)(12a_ + g2a+) , 
C = 2q;o(4q!_ +52(0!+ - 6)) , 

D = (^526 + 4a_ - i 520!+) (^a- + ^ 52"+) + 53ao(36 - 2a+) , 

and 

^i„t(x) = 2 ^ (p(a;^7.) + p(a;+)) a{a + 5„ ) . 

The restriction of this spin model to the polynomial space 5,„ (see Remark 4) yields the scalar 
elliptic CS model in jlTt provided and Q!_ = — 52a+/12. 
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Remark 5. Case 1 with a+ = yields the rational Calogero An spin model. Case 5 with a+a_ = 
is the model studied by Inozemtsev |^^, while for a+ = a_ = the hyperbolic Sutherland A]\f 
spin model is obtained. Case 6 with a+ — a- = = is the trigonometric Sutherland An spin 
model. The remaining potentials are new. 

In Cases 1-5, the potential is ES if a+ — 0. In Case 5, the potential is also ES for a_ = 0. The 
only ES potential in Case 6 is the trigonometric Sutherland potential (a+ — a_ = ao = 0). The 
remaining potentials, including all the elliptic potentials in Cases 7-9, are QES. 



Remark 6. In order to qualify as a physical wavefunction, the spin functions V'(x) constructed 



in Section III must vanish at all points in which the corresponding potential F(x) is singular. 



In addition, in Cases 1,2,3,5 the function ij) is required to be square-integrable over a suitable 
domain of K^. Both requirements impose certain constraints on the parameters and a defining 
the potential. A detailed analysis of the necessary and sufficient conditions on these parameters 
lies beyond the scope of this paper; see |Q for a complete solution of the analogous one-particle 
problem. However, it is not difficult in each case to provide sufficient conditions for the above 
requirements to hold. For example, in Case 5 the spin functions ?/; vanish at the singularities of the 
potential if and only if a > 0, and are square-integrable over the domain {x G : a;i > • • • > x^} 
provided < and a_ > 0. 



Remark 7. In Cases 1-6, the gauge factor is of the form 

n{^)^Y{[f{x-,)g{xt,)YX{h{x,), 



(34) 



i.<3 



with g — 1 or g = f. On the other hand, in the elliptic Cases 7-9 the gauge factor does not 
factorize as (34). This is consistent with the analogous result proved by Calogero for elliptic 
potentials in the scalar case m. 



V. CONCLUSIONS 



We have developed in this paper a systematic method for constructing new families of exactly 
and quasi-exactly solvable spin Calogero-Sutherland models. The key idea consists in relating the 
physical spin Hamiltonian to a general quadratic combination involving the usual A^v-type Dunkl 
operators (^, (^, and the new family of Dunkl operators introduced in Section || (Eq. (^)). 
Our approach goes beyond the Lie algebraic method extensively used in the scalar case, since the 
Dunkl operators (^ do not span a Lie algebra. However, they are invariant under the projective 
action of GL(2,R), a fact that is exploited in Section [V to classify all the potentials obtained by 
this method up to translations. 

The potentials constructed in this paper are all invariant under the An group consisting of 
permutations of the particles' coordinates and spins. A remarkable feature of the potentials in 
Cases 2-4 and 7-9 is their additional invariance under a change of sign of the spatial coordinate 
of any particle. Therefore, although these potentials are not invariant under the full Bn group 
of permutations and sign reversal of the particles' coordinates and spins, they are invariant under 
the restriction of the action of this group to the spatial coordinates. These models thus occupy an 
intermediate position between the usual spin CS models of An type and the fully _BAr-invariant 
rational and trigonometric spin CS models introduced by Yamamoto | |4^ . In fact, while Dunkl 
has recently proved the exact solvability of the rational Yamamoto model , there are no exact 
results for the eigenfunctions of its trigonometric counterpart. It is to be expected that a suitable 
extension of the method developed in this paper to the Bn case will yield new families of (Q)ES 
spin CS models of Bn type, including the trigonometric Yamamoto model. 
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VI. APPENDIX 

The formulae for the sums of the squares and the anticommutators of the Dunkl operators (| 
are given by 



E(^°)' - E (^'^' + (2 - + 2aE ^ (^.^a., - z]d.;) 



^' + ^? , ... ... a^^'., .. .. . W ^^^^ 



i j i<j ' ^ 

- « E (^^X {l-K,,)+aJ2{z, + z.fil - K,,) (A3) 



(A4) 



— 2am ZiZj + m(m — 1) , 

i<j i 

I Ei-^"' = E i'^^l + J(2 - b) a,,) + 2aE (^.a,, - z^a.J 

i i 'i'<j 

I E{^^' -^n - E i'f^. + (2 " -^s^J + 2aE ^ (^'^- - ^J^^.) 

i i i<ij 

^ «E (1 - i^.) - ^E'd - ^^.^.-) (^5) 

-— (l + a(Ar-l)), 

I E{^"' = E (-'^'. + 1(2 - ^) + 2«E 7^ (^'^- - ^1^^.) 

""E - + aE(^' + - ^^J-) (A6) 

m(2m — b) v-^ 

+ ^^E-^' 

where 

6 = 1 + m + a(Af - 1). 
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Here the symbol ^ ' stands for summation in i,j,k with i ^ j ^ k ^ i. The following identities 
are needed for the computation of the potential: 



E' ,.._..l_„., =>(^-i)(^-^). (M) 











Zi 






-Zj){z, 


- Zk) 












- Zk) 












- Zk) 









E r.. - ..L - .A = iN-m-2)Y: zu (Aio) 



{N -2)Y,{zi + Zi)\ (All) 



~ ^^^^^^ ~ ^*=) i<j 

5:-^ = liV(Ar-i), (A13) 

E74t = (^-1)E^- (A14) 

z^ 



E 7^ = - 1) E^i + E^^^^- (A15) 
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